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Abstract 



A search for the standard model Higgs boson in the H — > ZZ — > 2i2v decay channel, 
where £ = e or }i, in pp collisions at a center-of-mass energy of 7 TeV is presented. 
The data were collected at the LHC, with the CMS detector, and correspond to an 
integrated luminosity of 4.6fb^\ No significant excess is observed above the back- 
ground expectation, and upper limits are set on the Higgs boson production cross 
section. The presence of the standard model Higgs boson with a mass in the 270- 
440 GeV range is excluded at 95% confidence level. 
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1 Introduction 

The standard model (SM) of particle physics JlHll accommodates essentially all relevant exper- 
imental data. One of the remaining questions is the origin of mass for fundamental particles. 
Within the SM, vector boson masses arise from the spontaneous breaking of electroweak sym- 
metry II4U9J. The existence of the associated field quantum, the Higgs boson, has yet to be 
experimentally established. The discovery or exclusion of the SM Higgs boson is one of the 
main goals of the physics programme at the CERN Large Hadron Collider (LHC). 

To date, experimental searches for the SM Higgs boson have yielded null results. Limits at 95% 
confidence level (CL) on its mass {niu) have been placed by experiments at the Large Electron- 
Positron Collider (LEP), ma > 114.4GeV [lOJ, the Tevatron, ma i (162-166) GeV [llj, and 
ATLAS, mn i (145-206), {1\^-Tl^), and (340-450) GeV L12-14J. The primary production mech- 
anism for the Higgs boson at the LHC is through gluon fusion |[T5] - |26I with a small contribution 
from vector boson fusion (VBF) t27f|29l| . 

A search for the SM Higgs boson is presented in the H — > ZZ — > 2.17.V channel (where I refers 
to either e or yL), which is especially sensitive in the high-mass range 250-600 GeV. Results are 
reported from a data sample corresponding to an integrated luminosity of 4.6 fb^^ recorded in 
2011 by the Compact Muon Solenoid (CMS) experiment at = 7TeV. 

2 CMS Detector and Simulations 

A detailed description of the CMS detector can be found in Ref. [30 J. The key components of 
the detector include a silicon pixel and a silicon strip tracker, embedded in a 3.8 T solenoidal 
magnetic field, used to measure the momentum of charged particles. The silicon pixel and 
strip tracking system covers the pseudorapidity range |?/| < 2.5, where ?/ = — In [tan (0/2)], 
and Q is the polar angle of the trajectory of the particle with respect to the beam direction. It 
is surrounded by a crystal electromagnetic calorimeter (ECAL) and a brass-scintillator hadron 
calorimeter (HCAL). The ECAL and HCAL extend to a pseudorapidity range of |;/| < 3.0. 
A steel/ quartz-fiber Cherenkov forward detector (HE) extends the calorimetric coverage to 
\r]\ < 5.2. The calorimeters are surrounded by the muon system, used to identify muons and 
measure their momentum. The muon system consists of gas detectors placed in the steel return 
yoke of the magnet. 

The largest background to the SM Higgs boson signal consists of events in which a Z boson is 
produced in association with jets (Z -|-jets). The Z +jets cross section is five orders of magnitude 
larger than the expected production cross section for the signal. The other major backgrounds 
are top-quark production (tt — ?► 2£2i/2b and tW — > 2£2vh, and the diboson production (WZ — > 
3ev, ZZ leiv, and WW 2£2v). 

Several Monte Carlo event generators are used to simulate the signal and background pro- 
cesses. The H — > ZZ — ?► 2£2v signal and top-quark background events are generated by 
using the next-to-leading order (NLO) program POWHEG 2.0 1311 . The Z -|- jets and diboson 
backgrounds are simulated by using the MadGraph 5.1.3 generator |!32j. The diboson back- 
grounds are also simulated using the PYTHIA 6.4.22 generator for evaluating certain systematic 
uncertainties. Eor events generated using POWHEG and MadGraph generators, parton show- 
ering is simulated by using pythia with the Z2 tune which differs from the Zl tune described 
in Ref. [33J as it uses the CTEQ6 [34J parametrization for the parton distribution functions 
instead of the CTEQ5 [35| parametrization. The signal events are reweighted so that the trans- 
verse momentum (pj) distribution of the Higgs boson agrees with the next-to-next-to-leading 
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3 Event Selection 



order (NNLO) arid next-to-next-to-leading log (NNLL) prediction I|36ll37ll . The total cross sec- 
tion is taken from Ref. [38 J and is scaled by the H -> ZZ -> 2i2v branching ratio l|39H44| . 
The parton distribution functions (PDF) are modeled through the CTEQ6L [34J parametriza- 
tion at leading order and the CTIO parametrization |45| at NLO. The NLO contribution to the 
qq — > ZZ process is taken into account by reweighting the transverse momentum of the visible 
Z boson to match the prediction from the MCFM 6.0 program ||46| . A correction of 12% of the 
leading order qq ZZ cross section is included to account for the gg — ?■ ZZ process [47] . The 
detector response to the simulated events is modeled with GEANT4 |[30[[l8l and reconstruction 
and analysis are performed by using the same software used for data. 

3 Event Selection 

For Higgs boson masses considered in this analysis, the Z bosons from H — > ZZ decay are typi- 
cally produced with a substantial pj. Events are therefore selected to have two well-identified, 
isolated, opposite-charge leptons of same flavour (e+e^ or ^^}^^) with pj > 20GeV that have 
an invariant mass within 30 GeV window centred on the Z mass. The pj of the dilepton system 
is required to be greater than 55 GeV. In the electron channel, these events are collected by us- 
ing dielectron triggers, with thresholds of pj > 17 GeV and pj > 8 GeV for the leading and the 
other electron, respectively. The muon channel relies on a combination of single- and double- 
muon triggers. As instantaneous luminosity increased, the thresholds on the double-muon 
triggers changed from a requirement of pj > 7 GeV for each of the two muons to px > 17 GeV 
and Pt > 8 GeV on the leading and the other muon, respectively. The threshold for the single- 
muon trigger increased from pj > 17 GeV to pj > 24 GeV. The trigger efficiency for signal, 
for events selected through the full set of offline requirements is measured by using Z decays 
in data, and ranges from 95% to 97% in the muon channel and exceeds 99% in the electron 
channel. 

Muon candidates are reconstructed by using two algorithms, one in which tracks in the silicon 
tracker are matched to energy deposits in the muon detectors and another in which a combined 
fit is performed to signals in both the silicon tracker and the muon system [49 1. The muon 
candidates for analysis are required to be successfully reconstructed through both algorithms. 
Other identification criteria based on the number of measurements in the tracker and in the 
muon system, the fit quality of the muon track, and its consistency with the origin from the 
primary vertex are also imposed on the muon candidates to reduce the misidentification rate. 

Electron reconstruction also involves two algorithms [50], one in which energy clusters in the 
ECAL are matched to signals in the silicon tracker and another in which tracks in the silicon 
tracker are matched to ECAL clusters. The electron candidates used in the analysis can be re- 
constructed by either algorithm. More identification criteria based on the distribution of the 
shower in the ECAL, a matching of the trajectory of an electron track with the cluster in the 
ECAL, and consistency with origin of the track from the primary vertex are imposed on the 
electron candidates to reduce the misidentification rate. Electron candidates with an ECAL 
cluster in the transition region between ECAL barrel and endcap (1.4442 < < 1.566) are 
rejected. Additional requirements are imposed to remove electrons produced in photon con- 
versions in the detector material. 

Leptons produced in the decay of Z bosons are expected to be isolated from hadronic activity 
in the event. The sum of scalar transverse momentum depositions in the calorimeters and 
the transverse momenta of tracks in a cone of radius 0.3 in rj-cp space around each lepton, 
where (p is the azimuthal angle, is corrected by the contribution from the lepton and the ratio 
of this corrected sum divided by the lepton pj is required to be smaller than 15% (10%) for 
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muons (electrons). To correct for the contribution to the isolation sum from pile-up interactions 
(overlapping minimum-bias events from other concurrent proton-proton collisions), a median 
energy density (p) is determined event by event [51J. Then the pile-up contribution to the 
isolation sum is estimated as the product of p and the area of the cone in which the isolation 
sum is computed, and it is subtracted from the isolation sum to make it largely insensitive to 
pile-up. The combined reconstruction, identification and isolation efficiency is measured in 
data by using Z decays and ranges between 90% and 97% for muons, and between 70% and 
90% for electrons, depending on the pj and f] of the leptons. 

The high instantaneous luminosity delivered by the LHC provides an average of about 10 pile- 
up interactions per bunch crossing, leading to events with several possible primary vertices. 
The vertex with largest value of Pj for the associated tracks is chosen to be the reference ver- 
tex. According to simulation, this requirement provides the correct assignment for the primary 
vertex in more than 99% of both signal and background events. 

The presence of a large imbalance in transverse momentum in an event (£^^®^) is a fundamen- 
tal feature of the signal. The value of E^^^^ is the modulus of the E^'*® vector computed as 
the negative of the vector sum of the transverse momenta of all reconstructed objects identi- 
fied through the particle-flow algorithm, which aims to reconstruct all particles produced in 
a collision event by combining information from all sub-detectors [52 J. A large E^^^^ thresh- 
old is imposed to suppress the bulk of the Z+jets background, which contains little genuine 
£miss jj-^g region of large £^'®® is populated by Z+jets events in which the E^'^^ is largely 
due to jet mismeasurement. To suppress the background with £^^®^ arising from mismeasure- 
ment of jets, events are removed if the angle in the azimuthal plane between the E^^^^ and the 
closest jet with transverse energy Ej > 30 GeV is smaller than 0.5 radians. For events having 
no jets with Ej > 30 GeV, this requirement is imposed between E^^ and the closest jet with 
£t > 15 GeV. Jets are reconstructed from particle-flow candidates ||52l|53| by using the anti-kj 
clustering algorithm [54J with a distance parameter R of 0.5, as implemented in the FASTJET 
package [55, 56 J. 

Top-quark decays are characterized by the presence of jets originating from b quarks (bjets), 
which are tagged on the basis of impact parameters of tracks in a jet, relative to the primary 
vertex I|57ll58| . The top-quark background is suppressed by applying a veto on events having 
a b tagged jet with transverse energy greater than 30 GeV that lies within the tracker volume 
(|j/| < 2.4). To further suppress the top-quark background, a veto is applied on events con- 
taining a "soft muon" with pj > 3 GeV, which is typically produced in the leptonic decay of 
a b quark. The soft-muon veto along with the b-jet veto reduces the top-quark background by 
a factor of six. To reduce the WZ background in which both bosons decay leptonically, any 
event with a third lepton (e or ji) with pj > 10 GeV and passing the identification and isolation 
requirements is rejected. 



The search for the SM Higgs boson is performed by using a transverse mass (Mj) variable as 
the final discriminant in searching for an excess of events from the presence of the signal. The 
transverse mass is defined as follows: 



4 Analysis Strategy 
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5 Background Estimation 



where Pt{(1) and M{ii) are the transverse momentum and invariant mass of dilepton system, 
respectively. 

Two approaches are considered: a "cut-based" and a "shape-based" analysis. The same set 
of event selection criteria is used by both the analyses, except for requirements on the E^^^ 
and Mt variables. In both cases, requirements on £^^^ and Mj vary with the Higgs mass 
hypothesis, but are typically more relaxed for the shape-based analysis. 

In the cut-based analysis, the Higgs boson search is performed by looking for an excess of 
events above the standard model background expectation after applying selection criteria that 
are tuned for a given Higgs boson mass hypothesis. The £^'*^ and Mj selection listed in Table nl 
was optimized by using a genetic algorithm for rectangular cuts optimization, GARCON [59T7 
with the expected 95% CL limit on the cross section normalized to the standard model Higgs 
boson cross section as a figure of merit. Optimization was performed on six wjh points, and 
a smooth-fit interpolation for the cut values was performed to provide optimal cut values for 
every explored Higgs mass point. 

In the shape-based analysis, the Higgs boson search is performed by using a binned likelihood 
fit to the Mt distribution obtained after applying the selection given in Table [l] The range of Mj 
is varied as a function of nin to ensure that the expected signal distribution is fully contained. 
The shape-based analysis is found to have the best expected signal sensitivity and is considered 
as the main result in the paper. 

Table 1: Higgs boson mass-dependent selection for E^^^^ and Mj variables in the cut-based and 
shape-based analyses. 



Wh (GeV) 


250 


300 350 400 


500 


600 


Cut-based analysis selection 




> 70 


> 79 > 95 > 115 


> 150 


> 161 


Mr (GeV) 


[222, 272] 


[264, 331] [298, 393] [327, 460] 


[382, 605] 


[452, 767] 


Shape-based analysis selection 


£^'«« (GeV) 


> 70 


> 80 > 80 > 80 


> 80 


> 80 


Mr (GeV) 


[180, 300] 


[250, 350] [250, 400] [250, 450] 


[250, 600] 


[250, 750] 



5 Background Estimation 

For both types of analyses, the ZZ and WZ backgrounds are modeled using Monte Carlo simu- 
lation, and are normalized to their respective NLO cross sections. The remaining backgrounds 
are estimated using control samples in data. 

The Z+jets background is simulated from a control sample of events with a single photon pro- 
duced in association with jets (7 + jets). This obviates the need for using less reliable Monte 
Carlo simulation of the £J^'®® distribution arising from mismeasurement of jets. The 7 + jets 
cross section is much larger than Z + jets. The jets must have a minimum Ej of 15 GeV to re- 
duce contamination from processes that have a photon produced in association with genuine 
£^^*^, such as W{iv) + 7, W(£i/)+jets, where the jet is mismeasured as a photon, and Z{vv) + 7 
events. The kinematics and overall normalization of 7 + jets events are matched to Z + jets in 
data through an event-by-event reweighting as a function of the boson pj, and the number of 
jets with £t > 30 GeV, to account for the dependence of the £^'®® on the associated hadronic 
activity. Residual differences from pile-up are taken into account by reweighting events ac- 



5 



cording to the number of reconstructed vertices. This procedure yields an accurate model of 
the £^'®^ distribution in Z+jets events, as shown in Fig. [ij which compares the £^'®® distribu- 
tion of the reweighted 7+jets events along with other backgrounds to the E^^^^ distribution of 
the dilepton events in data. To compute the Mj for each 7+jets event, the value of pj{£i) is de- 
fined as the photon pj and the value of M{i£) term is chosen according to a probability density 
function constructed from the measured dilepton invariant mass distribution in Z+jets events. 




E^^' (GeV) 

Figure 1: The £^^^ distribution in data compared to the background simulation. The dielectron 
and dimuon channels are combined. Contributions from ZZ, WZ, non-resonant background 
and Z+jets background are stacked on top of each other. The E^'^® distribution in signal events 
for ?«H = 300 GeV is also shown. 

The background processes that do not involve a Z resonance (non-resonant background) are 
estimated by using a control sample of events with dileptons of different flavor (e^fi^) that 
pass the full analysis selection. This background consists mainly of leptonic W decays in tt, 
tW decays and WW events. Small contributions from single top-quark events produced from 
s-channel and i-channel processes, W+jets events in which the W boson decays leptonically 
and a jet is mismeasured as a lepton, and Z — > tt events in which t leptons produce light 
leptons and £^'^* are included in this estimate of the non-resonant background. This method 
cannot distinguish between the non-resonant background and a possible contribution from 
H — ?► WW — ?► 2£2v events, which are treated as part of the non-resonant background estimate. 

The non-resonant background in the e^e^ and final states is estimated by applying a 

scale factor (a) to the selected e^ji^ events: 

= OC^i X KSeyi, Nee = X Nep. 

This oi factor is computed from the sidebands (SB) of the Z peak (40 < M{££) < 70 GeV and 
110 < M{ie) <200 GeV) by using the following relations: 

_Nf _ Ag 

~ MSB' *e - ^sB' 
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6 Systematic Uncertainties 



where Nf^, N^^, and N^^ are the number of events in the Z sidebands in a top-enriched sample 
of e+e^, V^H ' ^rid e^fi^ final states, respectively. Such samples are selected by requiring £^^^ 
>70GeV and a b-tagged jet in the events. The measured values of a with the corresponding 
statistical uncertainties are = 0.58 ± 0.02 and ocg = 0.42 ± 0.02. 

Since the Mj distribution obtained for the non-resonant background is based on a small sam- 
ple of e^}i^ events scaled by a, this distribution is smoothed l,60J for use in the shape-based 
analysis. 

6 Systematic Uncertainties 

Systematic uncertainties include experimental uncertainties on the selection and measurement 
of the reconstructed objects, theoretical uncertainties on the signal and background processes 
which are derived from Monte Carlo simulation, and uncertainties on backgrounds determined 
from control samples in data. These are summarised in Table |2j 

The theoretical uncertainties on the predicted signal are evaluated by changing the parton dis- 
tribution functions and the QCD renormalization and factorization scales [38. 45. 61-64]. The 
theoretical uncertainties on the NLO cross sections of the ZZ and WZ backgrounds are es- 
timated by the same procedure of changing the parton distribution functions and the QCD 
renormalization and factorization scales by using the program MCFM 6.0. 

The cross section for on-shell Higgs boson production and decay is calculated in the zero-width 
approximation. Recent studies show that current Monte Carlo simulation does not describe 
the correct Higgs-boson line shape for nin > 300 GeV. These effects are estimated to corre- 
spond to an additional uncertainty of 10-30% on the theoretical cross section for 400 < mn < 
600GeVt38H651|66|- 

Various factors contribute to the experimental uncertainties that apply to processes derived 
from Monte Carlo simulation. These include uncertainties on the trigger efficiency and lepton 
selection efficiencies. The effect of lepton momentum scale and jet energy scale is also taken 
into account and is propagated to the evaluation of E^'^®. The uncertainty on the b jet veto is 
estimated by measuring the b-tagging efficiency in data from fully leptonic tt decays [58 j. The 
uncertainty due to the modeling of pile-up is evaluated by shiftng the mean of the distribution 
of the number of pile-up interactions up and down by one. 

The uncertainty on the estimate of the non-resonant background is dominated by the statis- 
tical uncertainty of the e^ji^ control sample. The uncertainty on the prediction of the Z+jets 
background is affected by any residual contamination of the 7+jets control sample from pro- 
cesses involving a photon and genuine £^'**. Consequently, this method provides only an 
upper bound on the prediction. To account for this contamination, the Z+jets background pre- 
diction is changed from zero to the estimate obtained from the 7+jets control sample in the 
limit-setting procedure. This is achieved by taking half of the estimate from the 7+jets control 
sample as the central value for the Z+jets background and assigning a 100% systematic uncer- 
tainty to this value. Since Z+jets contributes to only 10-15% of the total background, the large 
uncertainty on this background does not have a significant impact on the overall performance 
of the analysis. 

In the shape-based analysis, the difference in the Mj distribution with respect to an alternate 
Monte Carlo generator (madgraph versus PYTHIA) is used as an uncertainty in the modeling 
of the ZZ, WZ backgrounds. The uncertainty on the Mj distributions arising from the limited 
number of Monte Carlo simulated events is taken into account by scaling all the bins up and 



7 



down by the respective statistical uncertainties. This is done in a correlated manner to ob- 
tain two bounding distributions, within which the nominal distribution is changed. For the 
non-resonant background, where the nominal Mj distribution is obtained by smoothing the 
Mj distribution of e^ji^ events, the bounding distributions are taken from simulation to ac- 
count for any possible bias introduced by smoothing. The uncertainty for signal is assigned 
by reweighting the pj spectrum of the Higgs boson to the one obtained from changing the 
renormalization and factorization scales in the NNLO+NNLL calculation Il36ll37l . 

Table 2: Summary of systematic uncertainties on event yields of signal and background pro- 
cesses. Ranges refer to uncertainties that depend on the hypothesised value of jwh- 



Source 


Uncertainty [%] 


T 1 Tmi'noGitTT' 

L/Uil 111 HJJlL y 




PDF for Higgs production (gluon-fusion) 


8.2-10.4 


PDF for Higgs production (VBF) 


3.8-7.6 


PDF for qq ZZ 


4.8 


PDF for qq WZ 


5.5 


QCD scale for Higgs production (gluon-fusion) 


7.6-11.1 


QCD scale for Higgs production (VBF) 


0.4-2 


QCD scale for gg — > ZZ 


20 


QCD scale for qq ZZ 


6.2 


QCD scale for qq WZ 


8.5 


Higgs boson line shape 


10-30 


Trigger 


1 (for ee), 2 (for jiji) 


Lepton identification and isolation 


2 


Lepton momentum scale 


5 (for ee), 2 (for fifi) 


Jet energy scale 


1-1.5 


bjet Veto 


1-1.2 


Pile-up 


1-3 


Non-resonant background 


15-100 


Z+jets 


100 



-1 

for the 4.6 fb dataset are listed in Table 
based on the selections described in Table 



7 Results 

The event yields for the cut-based analysis, as a function of mn, and the estimated backgrounds 

3 The Mt distributions for the shape-based analysis 
1 for mii = 300 and 400 GeV hypotheses are shown in 
Fig.|2j No significant excess of events is observed over the expectation from the SM background 
and limits are set on the production cross section of the standard model Higgs boson as a 
function of Wh- 

The median expected and observed 95% CL upper limits on the cross section a x BR(H — > 
ZZ — > 2£2v) for 250 < nin < 600 GeV by using the cut- and the shape-based approaches are 
shown in Fig. |3] The measured ratio R of the 95% CL upper limit cross section a to the SM 
Higgs boson production cross section CTgM as a function of nin is shown in Fig. |4] These results 
are obtained by using the CLg method |l67H69l. For the cut-based analysis the SM Higgs boson is 
excluded in the mass range 310^65 GeV at 95% confidence level, while the expected exclusion 
limit in the background-only hypothesis is 305^70 GeV. For the shape-based analysis the SM 
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8 Summary 



Table 3: Background estimates, signal predictions, and observed number of events for an inte- 
grated luminosity of 4.6 fb^^ after the cut-based selection. The uncertainties include the statisti- 
cal and systematic components (in that order). The results for signal contain only the systematic 
uncertainties as the statistical uncertainties are negligible. For the non-resonant background, 
the 68% upper limit is quoted in the case where the background prediction is zero. 



(GeV) 


ZZ 


WZ 


Top/WW/ 
W+jets/Z^ TT 


Z+Jets 


Total 
Background 


Expected 
Signal 


Data 


250 


36.0 ±0.2 ±2.6 


24.0 ± 0.3 ± 2.0 


65.0 ±3.8 ±5.8 


15.0 ±15.0 


140.0 ± 3.8 ± 16.0 


22.0 ± 2.2 


142 


300 


23.0 ±0.2 ±1.7 


13.0 ±0.2 ±1.1 


18.0 ±1.1 ±3.0 


6.3 ± 6.3 


60.0 ± 1.1 ± 7.3 


21.0 ±2.1 


64 


350 


16.0 ±0.1 ±1.1 


7.0 ±0.2 ±0.6 


2.0 ±0.1 ±1.0 


4.1 ±4.1 


29.0 ± 0.3 ± 4.4 


21.0 ±2.5 


26 


400 


12.0 ±0.1 ±0.9 


4.6 ±0.1 ±0.4 


< 1.1 


2.7 ±2.7 


19.0 ±0.2 ±2.9 


17.0 ±2.0 


18 


500 


7.5 ±0.1 ±0.5 


2.0 ±0.1 ±0.2 


< 1.1 


1.4 ±1.4 


11.0 ±0.1 ±1.5 


7.4 ±1.3 


14 


600 


3.9 ±0.1 ±0.3 


0.8 ±0.1 ±0.1 


< 1.1 


0.6 ±0.6 


5.3 ±0.1 ±0.7 


2.9 ± 0.7 
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Figure 2: The Mj distribution for events passing mn selections for 300 GeV (left) and 400 GeV 
(right). The dielectron and dimuon channels are combined. 

Higgs boson is excluded in the mass range 270^40 GeV at 95% confidence level while the 
expected exclusion limit for the background-only hypothesis is 290-490 GeV. 



8 Summary 

A search for the standard model Higgs boson has been performed in the decay charmel H — > 
ZZ — > 2i2v in pp collisions at a/s = 7 TeV, using a data sample corresponding to an integrated 
luminosity of 4.6 fb^^. No significant excess is found above the background expectation. The 
presence of the SM Higgs boson is excluded for 270 GeV < wh < 440 GeV at 95% CL. 
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Figure 3: The median expected and observed 95% CL upper limits on the cross section a x 
BR{H — )■ ZZ — )■ 2£2v) for the Higgs boson masses in the range 250 — 600 GeV for the cut-based 
(left) and shape-based (right) analyses. 
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Figure 4: The ratio R of the 95% CL cross section upper limit c to the SM Higgs boson pro- 
duction cross section ctsm as a function of the Higgs boson mass /wh for the cut-based (left) and 
shape-based (right) analyses. 
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